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Abstract 

We have studied numerically the quantum-classical crossover in the escape- 
rate for an uniaxial spin system with an arbitrarily directed field. Using the 
simple quantum transition-state theory, we have obtained the boundary sep- 
arating the first- and the second-order crossover and the escape-rate in the 
presence of the transverse and longitudinal field. The results apply to the 
molecular nanomagnet, Mni2- 
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We consider an easy-axis ferromagnetic nanoparticle, or a molecular cluster, that has 
metastable or degenerate classical spin states. The direction of the magnetization may 
change due to two mechanisms. At sufficiently high temperature the rate of the magne- 
tization reversal T obeys the Arrhenius law, T ~ exp(— AU /ksT), with AU being the 
height of the energy barrier. At a temperature low enough to ignore the thermal activation, 
quantum tunneling comes into play with T ~ exp(—AU/hu) where uj is some temperature- 
independent frequency related to the shape of the metastable potential well. The crossover 
between thermal and quantum regimes has been intensively studied in nanospin systems 

This issue was first raised by Chudnovsky and Garanin, || who observed that the 
crossover in the spin Hamiltonian 7i = —DS^ — H X S X becomes sharp (first order) for 
h x {= H X /(2DS)) < 0.25 and smooth (second order) for 0.25 < h x < 1. In the exponential 
approximation, when only the transition exponent is concerned, the first- (second-) order 
crossover of the escpe-rate is characterized by the discontinuity (continuity) of dY{T)/dT at 
the crossover temperature, T . Subsequent calculations || rendered the boundary between 
the first- and the second-order crossover for the uniaxial model with a transverse and lon- 
gitudinal field. It was also pointed out ]T||2|^]] ^h.a.t ^ ne nonanalyticity of the rate for the 
first-order crossover disappears when quantum corrections to the exponential approximation 
are computed. The purpose of this Letter is to investigate the "sharpness" of the first-order 
crossover. This question is especially important in the light of a recent experimental evidence 
of the first-order crossover in Mn 12 |§. 

We shall focus on the crossover in the molecular magnet Mn 12 (S = 10) when the 
external magnetic field has both transverse and longitudinal components. In order to cal- 
culate the splitting of the excited states, we shall perform a numerical diagonaliztion of 
the Hamiltonian. Using numerical results, we then obtain the group of levels which make 
the dominant contribution to the thermally assisted tunneling. The full problem of the 
escape rate will be solved by mapping the spin problem onto a particle one ||. Summing 
contributions of all excited levels with account of quantum corrections, we will show that 
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the escape rate becomes analytic for the first-order crossover, but changes sharply around 
Tq, in contrast with the second-order crossover. We also obtain the boundary between the 
two types of the crossover and find that the first-order regime is greatly suppresed by the 
longitudinal field in accordance with experiment |5). 

According to the simple quantum transition-state theory, jltj the escape-rate in the 
temperature range T -C AU is given by 

r ( r)^ ? r m .e Xp (-^=), (i) 

where T m = u(E m )W(E m )/(27i), u)(E m ) is the frequency of oscillations at the energy E m , 
W(E m ) quantum transition probabilities, and Zq the partition function in the well. It is 
evident from Eq. (|l|) that the nonanalyticity of the rate is not expected around T for 
any escape process, even though the crossover can be sharp because of the exponential 
dependences of T m and thermal populations on the parameters. The rate for the first-order 
crossover, T/(T) is not analytic in the exponential approximation, when the escape rate is 
approximately given by the dominant term in the summation. However, the summation 
over all energy levels in Eq. ([!]) smoothens this non-analyticity. 

The model with an arbitrarily directed magnetic field is described by the Hamiltonian 

H = -DS 2 Z - H Z S Z - H X S X} (2) 

The zero field Hamiltonian has uniaxial symmetry with easy axis along z and hard plane, 
xy. H z is the longitudinal field which affects the height of the energy barrier, and H x is 
the transverse field which is responsible for quantum tunneling as well as for the reduction 
of the barrier. In the first approximation, this spin model describes the magnetic molecule 
Mn 12 . 

Within the thermally assisted model, the magnetization reversal occurs by quantum tun- 
neling from thermally excited magnetic levels at magnetic fields which bring into resonance 
the levels m and m! belonging to different potential wells. Denoting m to be the escape level 
from the metastable well, the resonance condition is that the levels m and m! have the same 
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energy when H z = kD where m! = —m — k and k is the bias index. The escape rate from 
any level is proportional to the product of the probability of the thermal occupation of that 
level and the probability of quantum tunneling from the level. In this respect, the dominat 
level for a given temperature is determined by the function Q 

«™> = 1^expp^), (3) 

where it is assumed that the sum of the linewidths of the m-th and m'-th levels substantially 
exceeds the level spacing U) m '{— -Ew+i — E m i). Here E rn {= —Dm 2 — H z m) is the energy 
level of the spin system without transverse field, and AE mm i is the splitting of the pair of in- 
resonance levels m and w! on the opposite side of the anisotropy barrier. It is seen from Eq. 
(H) that, since the escape rate decreases exponentially with decreasing temperature, larger 
longitudinal fields are necessary at lower temperature to produce an observable tunneling 
rate. 

Now, let us calculate the dominant level which maximizes f(m). In order to do 
that, we first need to find the range of m in the metastable well for a given transverse and 
longitudinal field, i.e., —S < m < m t < 0, when m t is the level which is near the top of the 
barrier and inside the metastable well. For H x = and H z = kD, simple analysis shows that 
m t = ~ — 1 where [x] gives the integer part of x. Since the height of barrier decreases 
with increasing H x , one expects m t < m°. To find the value of mt we express Eq. (|J) in the 
spherical coordinate and study the energy in the easy plane given by 

E(6,(j) = 0) = -DS 2 (cos 2 6 + 2h x sm0 + 2h z cos6), (4) 

where h x>z = H X>Z /(2DS). Writing the height of the barrier as AU = DS 2 (Au), m t is 
determined by the relation 
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m t = - 



2 +Sj(l-h z ) 2 -Au 



(5) 



where it is noted that, since An = (1 — h z ) 2 in the absence of the transverse field, we obtain 



m t = m". Also, m t = —S\Jh x {2 — h x ) at h z = 0. Numerical calculation of Am of Eq. 
leads to the results for m t (h x , h z ) shown in Fig. [[]. 



Next, we consider the frequency of the real-time oscillations, uj{E m ) at the energy E m . 
This quantity cannot be calculated with the use of the energy (Q) in the spherical coordinate 
system in which the physical quantity which is equivalent to the mass of the system in a one 
dimensional case is unknown. Thus, for the mapping of the spin problem onto a particle 
one, the corresponding energy-dependent frequency is of the form 

u (B) = a r|-Lr W -=^=)" 1 (6) 
'DJiW E - U(x) ' 



where X12 are turning points in the particle potential U(x) for a given energy E. Introducing 
the parameter p = (t/ sa( j — E m )/AU, the specific form of p at a small value of h x becomes, 

where p = at the top of the barrier. This gives m/S = —h z , i.e., m = —k/2 which is 
related to m t discussed previously. After some trivial manipulation of Eq. ([]), we obtain 
the dependence of uj on h x , h z , and m. As a result, the frequency of oscillation in f(m) can 
be numerically deduced by taking h x — > 0. In this limit the frequency also can be computed 
for the spin model having the energy levels E m = —Dm 2 — H z m as the inversed density 
of states, i.e., the energy difference between neighboring levels at energy E. This is simply 
given by u m ~ —D(2m + k) for S ^> 1. Now, in order to calculate the level splitting AE mm t, 
we first consider the formula of the perturbation theory || 

2DS m '~ m 



[{m — m — 



x 



(S + m')\(S -m)\ ll/2 



K'~ m - (8) 



(S -m')\(S + m)!_ 

This formula is compared with the results from the direct numerical diagonalization. As 
is illustrated in Table |, there is a disagreement between them for levels with m ^$ —6. 
However, noting that mt(h x — 0, h z — 0.1) = —2 is shifted to m t (h x = 0.05, h z = 0.1) = —6, 
the level m = —6 or —7 is important to study the type of the crossover and the escape-rate 
around To. Accordingly, we will perform the numerical diagonalization for the level splitting 
of f(m). 



Now, we are in a position to calculate the dominant level rrid which is determined by the 
maximal value of the function ([|) for a given temperature. Within the thermally assisted 
tunneling model the system tunnels through the levels between the bottom and the top. In 
this process maiT) behaves in two different ways. One way is that rrid changes continuously 
from —S to m t , whereas the other way is that it performs some discontinuous jump. We 
call the former - the second-order crossover and the latter - the first-order crossover. As 
is shown in Fig. |2], for the resonant field h z = 0.1 we have the first-order crossover for 
h x = 0.05,0.1 and the second-order crossover for h x = 0.15. Also, for h z = 0.4 ( H z ~ 3.28 
Tesla in Mni2 ), the abrupt shift occurs at h x = 0.04 and the corresponding dominant level 
changes by 2 (m^ = —10 and rrid = —8) in the range of temperature (~ 0.1 K - ~ 1 K ). 
Strikingly, this feature is observed in a recent experiment, in which the step positions 
shift abruptly at low temperature and high magnetic field. Employing these schemes in 
the whole range of h z , we obtain the phase boundary for the values of the transverse field, 
which is shown by the symbols in Fig. |3|. In the quasiclassical method, the order of the 
quantum-classical escape-rate crossover was determined by the sign of the coefficient in the 
expansion of the imaginary-time action near the top of the barrier ||. In the perturbation 
method, the behavior of the result (Q) is inserted into the calculation of the dominant level 
rrid in the metastable well, whose behavior determines the type of the order. Using this 
method, the first-order crossover is found to be suppresed as compared to the quasiclassical 
method. This was noticed in Ref. 0. However, as discussed previously, the corresponding 
level splitting is not quite correct, especially, near the top of the barrier, in which range 
the perturbation fails. The correct calculation based on the diagonalization method shows 
that the first-order regime is even more suppresed as compared to the perturbative results. 
For example, in the unbiased case, the phase boundary becomes at h x = 0.114, 0.139, and 
0.25 for the diagonalization, perturbation, and quasiclassical method, respectively. In Fig. 
|3], there is no data point beyond h z = 0.8. The reason is that rrit = —9 in this region 
and thereby it is meaningless to ask whether the shift is continuous or discontinuous in this 
region. 
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The values of the crossover temperature T at the phase boundary between first- and 
second-order crossover, which have been obtained by the diagonalization method described 
above, are shown in Fig. |j. Comparing this with two other methods, the result based on 
the diagonalization method is the smallest in the whole range of the bias field. For example, 
in the unbiased case, we have T^ c) /(DS) = 0.117, 0.124, and 0.137 for the diagonalization, 
perturbation, and quasiclassical method, respectively. 

Now, for the relaxation at resonance we will present the results of numerical calculation 
for the escape rate, T(T). As is shown in Fig. |5| (as, e.g., for h z = 0.1), the escape rate in 
the first-order region (h x = 0.05 or 0.1) changes sharply, as shows the comparison with the 
escape rate in the second-order one (h x = 0.15 or 0.2). Also, for h z = 0.4 we can clearly 
distinguish the behavior of the rates in two different regimes, e.g., h x = 0.04 and 0.1, and 
its trend continues in the whole range of the bias field. The origin of these behaviors is 
that, as described above, changes discontinuously around T for the first-order crossover, 
while it does continuously for the second-order one. In other words, since the rate ([I]) is 
sensitive to the change of m, the abrupt jump in the first-order crossover induces sharp 
increase, e.g., T(h x = 0.05, h z = 0.1)/T{h x = 0.2, h z = 0.1) ~ 10 6 at T/(DS) ~ 0.12, and 
T(h x = 0.04, h z = 0A)/T(h x = 0.1, h z = 0.4) ~ 10 5 at T/(DS) ~ 0.11. 

For the uniaxial spin model considered in this paper, both type of crossover can be real- 
ized in the molecular magnet Mn 12 , and the situation can be controlled by the longitudinal 
and transverse field. For Mni 2 , D ~ 0.55K, and AU = DS 2 ~ 55K at h x = h z = 



TT|,[T2[. The critical field in the x- or z-direction is H xc = H zc = 2DS / (que) — 8.2T, 



and the longitudinal field for the resonance is H = D/^qhb) — 0.41T. At these fields, the 
magnetic relaxation in Mni2 can occur on measurement time scales, and gives rise to the 
different behavior of the dominant level and the rate in magnetization depending on H x 
and H z . Furthermore, in the unbiased case the first-order crossover can be observed in the 
field range < H x < 0.93 Tesla and the crossover region occurs at the temperature range 
~ 0.1 K < T <~ 1 K. This field range decreases with increasing the bias field. Even though 
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we presented the results of only for h z = 0.1 and 0.4 in Fig. |^, we have found that the 
levels which, at a certain temperature, change by more than 1 are m = 7, 8 or 9 depending 
on the value of the resonant field. 

In conclusion, we have studied the quantum-classical crossover of the escape-rate of 
a uniaxial spin model with an arbitrarily directed field. Employing the diagonalization 
method, we have obtained the dominant level for the thermally assisted tunneling and 
dependence of the escape rate on temperature. In comparison to the previously studied 
models, the first-order region is greatly suppressd, but still observable in molecular magnets. 
It is also found that the first-order crossover is fairly sharp while the second-order the 
crossover is smooth. This is found to be strongly related with the discontinuous (continuous) 
jump of md in lower (higher) field h x . These results have been applied to the high-spin 
molecule, Mn 12 . They are also relevant to the study of nanoparticles. 
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TABLES 

TABLE I. Comparison of the splitting (||) with the one from the numerical diagonalization for 
S = 10, h x = 0.05, and h z = 0.1. 



m 


pert. 


diag. 


-6 


0.156 


8.53 x 10" 


-2 


-7 


1.01 x 10~ 3 


7.38 x 10" 


-i 


-8 


2.72 x HT 6 


2.26 x 10" 


-6 


-9 


3.14 x HT 9 


2.78 x 10' 


-9 


-10 


1.40 x 10~ 12 


1.28 x 10- 
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FIGURES 



FIG. 1. rrif/S vs. h x for a given value of h z . The range of h z is ~ 0.9. 

FIG. 2. m d vs. T(= T/{DS)) for fe z = 0.1, where h x = 0.05 (a), 0.1 (b), and 0.15 (c). Inset: 
h z = 0.4, where h x = 0.04 (a), 0.06 (b) and 0.1 (c). 

FIG. 3. Phase boundary between the first- and the second-order crossover obtained by the 



quasiclassical (b), the perturbative (c) and the diagonalization method (d). The critical field [13] 



iJiJc + hxi 3 = 1) is represented in (a). 

FIG. 4. Crossover temperature Tq (= Tq C ^ /(DS)) at the phase boundary between first- and 
second-order crossover, based on the quasiclassical (a), the perturbative (b) and the diagonalization 
method (c) 

FIG. 5. f (= r(T)/r(0)) vs. t(= T/(DS)) for h z = 0.1 , where h x = 0.05 (a), 0.1 (b), 0.15 
(c), and 0.2 (d). Inset: h z = 0.4, where h x = 0.04 (a), 0.06 (b) and 0.1 (c). 
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